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ABSTRACT 



Experimental tables of availabilities at time t are 
obtained for a device whose performance is described by an 
alternating renewal process with a finite number of failure- 
renewal cycles, corresponding to having a finite spares 
backlog. Failure and repair rates are assumed to be constant, 
and attention is restricted to cases in which the repair rate 
is larger than the failure rate. 
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I. INTRODUCTION 



The most commonly encountered working definition of the 
"availability" of a device, 

MTTF 

(1.1) Availability = 

MTTF + MTTR 

represents the long-term or steady-state probability that 
the device will be found in an "up" or functioning condition 
when two specific conditions are satisfied. One condition is 
that there is an alternation of failure and repair cycles 
in which times to failure and times to repair are independent 
realizations from some failure and repair distributions 
satisfying minimal regularity conditions. The second, and 
here most important condition, is that the alternation of 
failure and repair continues indefinitely, so that the 
performance of the device is described by a standard alternating 
renewal process. 

For many equipments, the second condition cited above 
implies access to an infinite backlog of spares. In many 
operational contexts this sort of spares support cannot be 
realistically assumed. 

If availability is considered to derive from an alter- 
nating renewal process with a finite number of cycles, 
corresponding to a finite backlog of spares, then expressions 
for availability become complex as compared with equation (1.1). 
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This thesis is devoted to computational experiments with 
some "finite spares" availability expressions. The end 
objective of such experiments is to be able to determine the 
circumstances in which equation (1.1) furnishes an adequate 
approximation, or alternately to be able to provide compu- 
tationally feasible alternatives to its use. 
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II. MATHEMATICAL MODEL 



The mathematical model on which the usual expressions 
for the availability of a device are based is an alternating 
renewal process with an infinite number of failure-repair 
cycles. In situations where repair requires replacement of 
a failured device by a spare, this corresponds to having an 
infinite number of spares. The model studied here is modi- 
fied to allow only a finite number of failure -repair cycles, 
corresponding to having a finite number of spares. 

The simplest assumptions about failure and repair times 
are made; failure rates are constant, and repair rates are 
constant. Only those processes that begin with a functioning 
device installed are considered. 

In greater detail, the failure-repair process considered 
is as shown in figure 2.1, 
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..Figure 2.1 Failure-repair process. 



where n is the number of spares, Tq is the time to failure 

for the original device, T^, T^, T n are the times to 

failure for the n spares, and S Q , S^, S n-1 are t * ie 

times to replace the original device and the first n-1 spares. 

It is assumed that T n , S n , T ,, S ,, T are inde- 

u u n-l n-1 n 

pendent random variables, and that T Q , , T n are 

exponentially distributed with failure rate A, while 

S Q , S i are exponentially distributed with repair 

rate n. 

The availability at time t, A (t) , of the original 
device, supported by its backlog of n spares, is the 
probability that the process shown in Figure 2.1 is in an 
"up" condition at time t; i.e., that at time t either the 
original device or one of its spares is installed and still 
functioning. 

The increment in availability at time t due to the 
th. 

k spare, I^Ct) is defined by 

(2.1) I k (t) * A k (t) - A k _ 1 (t) k = 1, , n. 

Before proceeding to a derivation of expressions for 
I n (t) and A n (t) in a general case, two special cases are 
considered; repair rate n equal to infinity, and repair 
rate n equal to failure rate A. These are boundary cases 
for the cases of likely practical interest, in which it is 
reasonable to expect that repair rate will exceed failure 
rate. 
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In any case, Ag(t) , availability at time t with no 
spares is given by 

(2.2) A Q (t) = P [Tq > t] = e' Xt , t > 0. 



In the following sections, it will be convenient to let 



(2.3) 



U = T, 
n ( 



+ T 



n’ 



V n = S 0 + 



+ S 



n’ 



W = U + V 
n n n 



•■(S 0 *T 0 J * * (S n ♦ T n ). 

A. REPAIR RATE EQUAL TO INFINITY 

The simplest case is the one in which no time is required 
to repair a failed unit, provided a spare unit is available. 
In this case the contribution of the first spare is 



(2.4) 



r z 

qtt) - Ptq > t-s (Uq 

Jo 



s] fy (s) ds 
o 



where 

f n (s) = Xe“ Xs , s > 0 , 
u 0 



is the gamma {1,A} density. Thus 



(2.5) 



Ij_(t) 




‘ e -X(t-s) 



\e 



-Xs 



ds 
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, -As 

= Ae t, 



and the availability of a system having one spare is 



( 2 . 6 ) 



AjCt) = A 0 (t) + I 1 Ct) 



= e' Xt (l + At) 



The contribution of the second spare is 



(2.7) I 2 (t) = I P[T 2 > t-s|U 1 = s]fy (s) ds , 



where 



,2 -As 

f u, = — rfij - > s i 0 » 



is the gamma {2, A} density. Thus 



( 2 . 8 ) 



I 2 (t) - 



A (t-s) A 2 se' Xs ^ 

6 - T(2) as 



- ( A t) e jj - 



and the availability of a system having two spares is 



(2.9) 



A,(t) = e" Xt (l + At + 



th 

Generally, the contribution of the n spare is 
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t 



( 2 . 10 ) 



! n U) - 



P [T > t-s I U 
L n 1 



n-1 



= s]f 



U 



n-1 



(s) ds 



where 



L U 



n-1 



,n n-1 - A.s 

^ ■ ^TTnT i 0 • 



is the gamma (n, A} density. Thus 



( 2 . 11 ) 



I n (t) = 



- (t-s) .n n-1 -As 
e X s e 



TW 



ds 



- (\t) n e' Xt ij, , 



and the availability of a system having n spares is 
(2.12) A n (t) = e' At (l + At + + + ) . 

Availabilities A (t) obtained from equation (2.12) are 
shown in table 1. In this table n represents the number of 
spares, and the contribution I (t) of the n l spare can be 
found by subtracting A ^(t) from A (t) . The availabilities 
shown in the last column are for an alternating renewal pro- 
cess with an infinite number of failure-repair cycles. This 
corresponds to having an infinite number of spares. 

B. REPAIR RATE EQUAL TO FAILURE RATE 

The failure-repair process considered is that in which 
failure times and repair times are exponentially distributed 
with equal rates, i.e., A = n . 
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In this case, the contribution of the first spare is 



(2.13) I^t) =J P[T 1 > t-s|W Q = s]f w ^(s) ds , 



where 



, 2 _ Xs 

( s ) = ~ T TT T > s 1 0 > 



tot 



is the gamma (2, A} density. Thus 



(2.14) 



I,(t) = 



a/-, \ \2 -As 

- A (t- s) A se 

e TXH dS 



l'(2) 2 ’ 



and the availability of a system having one spare is 



A x (t) = A 0 (t) + I-l Ct) 



= e 



" xt (i + -^r-) 



1” Vi 

Generally, the contribution of the n spare is 
(2.15) I n (t) = f P[T n > t-s|l» n . 1 = s]f w (s) ds , 

Jo n_1 



where 



,2n 2n-l -As 

f '“ W ' ~ 1‘ ( 2n) 6 • 5 i 0 • 



L w 



n- 1 



13 



is the gamma {2n,A} density. Thus 




t 



(2.16) I n (t) 




l 



t 



s 



2n- 1 



ds 



(2n-l) ! 




2n -At 



and the availability o£ a system having n spares is 



(2.17) A n (t) = e' Xt (1 + + 




C. REPAIR RATE GREATER THAN FAILURE RATE 

The failure-repair process considered in this section 
is visualized as one in which the repair rate is greater 
than the failure rate. This influences the format in which 
the results are displayed. 

til 

In this case, the contribution of the n tx spare is 



r z 



(2.18) I n (t) 



p [T n > t-s|W n _ 1 = s] f w (s) ds , 



n-1 



where 




f U (t-s) ds , 



n-1 
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Table 2 A (t) for the case nt 
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i 



and 



n 



£ v (s) = fJKJ sn " le " nS , s t 0 > 



n-1 



is the gamma {n,r|} density, 
while 



n 



f U = TJKY (t-s) 11 ' 1 e' X(:t ' s) , t-s > 0 , 



is the gamma {n,A} density. Thus 



(2.19) I n (t) = 



.■A (t-s) [ n n „n-l -nu 1 A n . .n-1 

rTnl e rw c uj 



-A(t-s) , , 

• e v J du ds • 



Inverting the order of integration, equation (2.19) 
becomes 

>t r t 



(2.20) I n (t) = 




-A(t-s) ri n ,,n-l -riu A n , .n-1 

e riHl" 6 HnT C5 ‘ U) 



r X ( 5 - u ) ds du 



= e 



-At n n A n 



UCnjr Jo 



u n-l e -(n-A)u 



-J 



r z 



u 



(s-u) n ^ ds du 



Let v = s - u. Then 



■J 



r t 

0 



(s-u) n_1 ds = 



p t-u 
0 



v 11 ' 1 d V = 



n 



n 
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Thus equation (2.20) reduces to 



n - 71 '! t n n X n / ,.n-l,. ,,^n^-(n-X)u 

(2.21) I n (t) - e u (t-u) e 

J 0 



and the availability of a system having n spares is 



(2.22) A ft) = A n (t) + I-i (t) + + I (t) 



n 



n 



Note than when n = A., equation (2.21) reduces to 
equation (2.16), since then 



T -Xt I 2 " 

I n (t) e TTn)r(n*l) 

J 



u n_1 (t-u) n du 



Let u = tv. Then 



r t 
0 



u n_1 (t-u) n du 



r 



(tv) n-1 {t(l-v) } n tdv 



. 2n / n-1,.. ,n , 

t v ( 1 -v ) dv 



_ „2n r(n)r(n+l) 
r T (2n+l) 



du , 
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so that 



e" Xt X 2n 
I n ( ' t) f (nj T (n+1 J 

. (U) 2 V U 



2n r(n)r(n+l) 

i' (2n+l J 
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III. APPROXIMATION OF THE MATHEMATICAL MODEL 



In section II, we derived a mathematical model for 
availability. 

In this section, we discuss methods of approximation to 
obtain numerical values of availability. 



A. EXPONENTIAL EXPANSION APPROXIMATION 

The integral in equation (2.21) can be approximated by 
expanding its exponential term, i.e., 



(3.1) 



(n-*) = 1 



(n-x)u * 



u 



Thus the integral becomes 



(3.2) [ u n ~ 1 ( t - u) n {l- (n-A.) u + u2 



u n ^(t-u) n du 



rt 



- (n-x) 



U n (t-d) n du 



x (n-A.) 

TT~ 



r t 



,„.n+l .n , 
u (t-u) du 



} du 



We know that , 
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u n * 1 (t-u) n du 



2n r (n) r (n+1) 
T(2n+lJ 



Thus equation (3.2) becomes 



(3.3) 






, e' At q n X n r „2n r(n)r(n+l) 
T (n) r (n+l) r(2n+l) 



- (q-x) 



2n+l r(n+l)r(n+l) 
n7n+7) 



+ (q-x) 2 t 2n+2 r (n+2) r (n+i) 
n i' (2n+3 J 



} 



= e' Xt (qt) n (At) n , T(n) 
Hn) l r(’2n+l) 



(qt-At) 



r(n+i) 

IT2nV2) 



_ (qt-At) 2 r (n+2) 

T! r ( 2n + 3J 

- } 

Computational experiments, with the approximation 
represented by equation (3.3) have indicated unsatisfactory 
convergence behavior when (qt-At) is large, a case of some 
practical interest, and so this approach was not pursued. 

B. SIMPSON’S RULE APPROXIMATION 

The integral in equation (2.21) can be approximated by 
using Simpson's rule. 
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Then 



Let v = :jr 



(3.4) 



n- 1 . ~(r|-X)u j 

u (t-u) e v yu du 



n-1 



w (t(i-v)}" e ' Cn ' x)tv t dv 



r 1 

= t 2n v^d-V)” e' (,lt '* t)V dv . 

J 0 



Thus equation (2.21) becomes 



(3.5) I n (t) 




(nt) n (Xt) n 
i’TnJ i 1 (n + l) 




y 71 (1“ v) 



n 



e - (nt-Xt) v 



dv . 



Now Simpson's rule can be applied to the integral in 
equation (3.5) to obtain numerical values of availability. 
Simpson’s rule as applied is 



(3.6) 



b 

f(x)dx - j(y 1 + 4y^2y 3+ 4y 4+ .... + 2y m , 1+ 4y m .y mtl ) 



where 




9 



and 



Xi = a, X 2 = a+h. 



X m + 1 



a+mh = b , 
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while 



«Vi> 



yj - fCXj), y 2 ■ ftx 2 ). 



’ ^m+1 
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IV. TABLE 



In this section experimental tables of availabilities 

A ft) are shown, which were obtained from the mathematical 
n 7 

model evaluated using Simpson's rule, i.e., using equation (3.6) 
to evaluate equation (3.5) with m = 500, = 0.0001 and 

X 501 =0 .9999. The reason for choosing m = 500 is that compu- 
tational experiments with choices of m greater than or equal 
to 500 gave a stable and accurate result. 

In these tables n represents the number of spares, 
and the contribution of the n l spare I n (t) can be found by 
subtracting A ^(t) from A n (t). 

The availabilities shown in this section, Table 3-12, 
are for cases in which nt > Xt. 

The availabilities shown in the last column are for an 
alternating renewal process with an infinite number of failure- 
repair cycles . This corresponds to having an infinite number 
of spares. 
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Table 8 A (t) for nt = 120 
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Table 9 A (t) for nt = 140 



8 



/■ — \ 


04 


h- 




Cn 


CTi 


o 


04 


LO 


cn 




4-> 


00 


vO 


H* 


04 


LO 


cn 


04 


LO 


CO 


04 


v — / 


CT> 


Cn 


cn 


CTi 


CO 






vO 


LO 


LO 


8 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


CTi 


cn 


cn 


< 























CTi 



CO 



vO 



to 



4-» 

v — / 

o 

rH 

< 






00 

< 



o- 

< 



vO 



to 



<P 

v — / 

p- 

< 



o 


rH 


H- 


o^ 


cn 


O 


CO 


cn 




cn 


00 


vO 


H- 


04 


LO 


cn 


o 




to 


o 


04 


cn 


cn 


cn 


CO 


o^ 




LO 


to 


cn 


cn 

• 


cn 

• 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


00 


o 


rH 






cn 


LO 


H* 


LO 


H* 


o 


CO 


vO 


H* 


04 


LO 


CO 




LO 


to 


LO 


cn 


cn 


cn 


cn 


CO 




vO 


H* 


o 


to 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


CO 


o 


rH 




o- 


CO 




o 


cn 






CO 


vO 




04 


LO 


vO 


CO 


o- 


o- 


o- 


cn 


cn 


cn 


cn 


OO 


o- 


LO 


rH 


H* 




cn 


cn 


cn 


cn 




cn 


Oi 


cn 


CO 


O'* 


o 


rH 


h- 


o^ 


04 


vO 


LO 


vO 


oo 


to 


CO 


vO 


H" 


04 


LO 


o 




to 


VO 


VO 


cn 


cn 


cn 


cn 


CO 




to 


vO 


LO 


04 


cn 


cn 


cn 


cn 


cn 


cn 


cn 


CO 


o- 


vO 


o 


rH 




vO 


LO 


rH 




O 


04 


cn 


CO 


vO 




04 


04 


to 


to 


O 


o- 


CO 


cn 


cn 


cn 


cn 


CO 


LO 


CO 


o- 


04 


o-* 


cn 


o> 


cn 


cn 




cn 


CO 


o- 


vO 




o 


rH 


to 


04 


H- 


o 




CO 


vO 


H" 


CO 


vO 


H- 


04 


04 


cn 


CO 


o 


o- 


H* 


cn 


cn 


cn 


cn 




o 


CO 


to 


vO 


04 


cn 




cn 


cn 


cn 






VO 




to 


o 


o 


LO 


LO 


rH 






o 


CO 


o-* 


CO 


vO 


to 




O 


LO 


CO 


vO 


rH 


o-* 


cn 


cn 


cn 


CO 




rH 


to 


LO 


O 


CO 


cn 


cn 


cn 


cn 


cn 


00 


VO 




to 


rH 


cn 


vO 


CO 


LO 


o 




LO 


CO 


to 


H* 




H" 


o- 


LO 


LO 


04 


to 


LO 


o 


CO 


cn 


cn 


oo 


O- 


LO 


LO 




o- 


vO 


CO 


CT> 


cn 


cn 


cn 


co 


vO 




04 


rH 


O 


o 


vO 


04 


CO 


LO 


vO 


o 


H- 


04 


vO 


vO 


04 


vO 


vO 


CO 


CO 




o> 


LO 


rH 


cn 


CO 


LO 


rH 




04 




04 


vO 


to 


cn 


cn 


cn 


cn 


vO 




04 


rH 


O 


O 


04 


to 


LO 


to 


vO 


rH 


rH 


H* 


cn 


vO 


04 


oo 


LO 


LO 




rH 


to 


rH 


o^ 




O- 


o 


04 


to 


o 


O 


cn 


H- 


rH 


O 


cn 


cn 


CO 




H- 


04 


o 


O 


O 


o 



4-> 

r-< 



LD 

04 



LO 



LO 



04 



to 



LO 



vO 



30 





r — \ 




a* 


to 


CO 


O' 


O 


vO 


o- 




i— i 




P 


CO 


so 


un 


to 


O* 


i— 1 


LO 


a\ 


to 


CO 


8 


v / 


o* 


o* 


cn 


cr> 


CO 


CO 


o* 


vO 


vO 


LO 


8 


a\ 


CT* 


cn 


cn 


a* 


a* 


o> 


C7* 


G* 


cn 




< 
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P 
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o 


vO 


i— i 
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Table 11 A (t) for nt = 180 
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V: SUMMARY AND CONCLUSIONS 



Certain computational approaches were tried for obtain- 
ing availabilities for a device supported by only a finite 
backlog of spares, using the simple assumptions that failure 
and repair rates are constant. Real failure and repair distri 
butions may be more complex, but the case considered is a 
good case for initial computational experiments. 

Of the two approaches tried, neither proved entirely 
satisfactory in obtaining availabilities in a way that is 
fast and suitable for use with small-scale computational 
facilities, e.g.. hand-held calculators. Also neither was 
effective, over the entire range of failure rate and repair 
rate combinations that might be of interest. 

Since an easily used, readily accessible , way to assess 
the impact of finite spares backlogs on availability is 
desirable in many mission planning contexts, further computa- 
tional approaches should be tried. 

The tables presented in section IV give availability 
values with which the result of such experiments can be 
compared. 
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